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The occurence of spin-polarization at Zr02, AI2O3 and MgO surfaces is proved by means of ab- 
initio calculations within the density functional theory. Large spin moments, as high as 1.56 hb, 
develop at O-ended polar terminations, transforming the non-magnetic insulator into a half-metal. 
The magnetic moments mainly reside in the surface oxygen atoms and their origin is related to the 
existence of 2p holes of well-defined spin polarization at the valence band of the ionic oxide. The 
direct relation between magnetization and local loss of donor charge makes possible to extend the 
■ magnetization mechanism beyond surface properties. 
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r/3 ' When dimensions are reduced to the nanoscale, we are faced to a new understanding of the physical properties of 
1 1 , matter: bulk insulators and semiconductors exhibit metallic surfaces [1], non-magnetic materials get spin polarization 
when forming nanoparticles [2], unstable bulk structures exist in ultrathin film form [3, 4], etc. At the origin of these 
phenomena are the reduced dimensionality and the enhanced role of the surfaces or boundaries in the final properties 
of the system. Together with its inherent fundamental interest, this has important technological consequences, aus- 
picating the birth of new technologies [5, 6]. Special attention is devoted to magnetic low dimensional structures, in 
particular as sources of spin current in the emerging field of spintronics [5] . 

In this letter, we report on the existence of large magnetic moments and half-metallicity at the O-rich surfaces of 
ceramic oxides, focusing on ZrC>2, AI2O3 and MgO. These are non-magnetic ionic insulators widely applied in bulk 
and thick film form, that have also been grown as ultrathin films and nanometric grains [7] . Their electronic structure 
' can be roughly described as a valence band formed by the filled O 2p orbitals and a conduction band formed by 
l— ~~ ', the empty metal levels. When a M x O y unit -M being the metal donor and x,y accounting for the particular metal 
to oxygen ratio- is broken to form the surface, the loss of coordination of the surface O atoms originates 2p holes 
, in the valence band of the oxide. Our results show that this generates high spin moments at the topmost O layer, 
which induce magnetization at the adjacent planes and, remarkably, alter the electronic structure of the oxide from 
r insulating to half-metallic. 

Very recently unexpected ferromagnetism has been measured in thin films of undoped non-magnetic oxides, like 
Hf02 and Zr02, possibly assigned to the presence of lattice defects concentrated at the film interface. [8, 9]. Here we 
■ prove the existence of a magnetization mechanism rooted in the loss of donor charge of the O atoms, something that 
can also occur in films with cation vacancies. This provides an explanation for the origin of the magnetic moments 
• of these so called 'd-zero' ferromagnets. Furthermore, we predict that this kind of magnetism can be extended to a 
wider class of non-magnetic oxides, like MgO and AI2O3, where the cation is not necessarily a d metal. 

The presence of defects has already been found to transform certain bulk ionic materials from insulating non- 
magnetic to ferromagnetic: either by adding transition metal dopants [10], as occurs with Ti02, Sn02 and ZnO, or 
by a kind of molecular Hund's rule coupling for the solid state [11], as in the case of dilute divalent Ca vacancies 
in CaO and dilute neutral 3$ vacancies in the hexaborides. [12]. Although the magnetization mechanism reported 
O ■ here has an electrostatic origin similar to the one proposed in bulk CaO, new ingredients are present: the localization 
of the spin moments close to the surface, and the extended character of the polarized 2p valence band states. This 
V '"j ■ suggests that intrinsically low-dimensional structures of the oxides considered, such as grains or ultrathin films, can 
rS I be magnetic. 

Our results are obtained from first-principles spin-polarized calculations within the DFT (density functional theory) 
under the GGA (generalized gradient approximation) for exchange and correlation [13]. For highly correlated systems 
such as transition metal oxides (TMO), this approach predicts too small magnetic moments and inaccurate magnetic 
energies due to its failure in describing the localized nature of the d electrons [14]. However, this is not the situation 
for the ceramic oxides considered here, with a valence band formed by O 2p orbitals. In fact, there are no valence d 
electrons in MgO and AI2O3, while it is widely recognized the success of the DFT-GGA formalism, with and without 
self-energy corrections [15], to determine the structural and electronic properties of Zr02, which are in excellent 
agreement with experiments [16]. We use the SIESTA package [17] with basis sets formed by double-zeta polarized 
localized numerical atomic orbitals. More details about the conditions of the calculations can be found elsewhere [18]. 
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We have studied binary insulating oxides of different crystal structures: MgO, AI2O3 and ZrC>2, for which few 
calculations have considered the possibility of spin polarization at clean surfaces or in defective structures [4]. MgO 
adopts the rock-salt lattice, the corresponding (111) surfaces consisting of a stacking of alternated pure Mg and O 
planes. CK-AI2O3 can be viewed as an hexagonal close-packed array of O atoms with the Al atoms occupying two 
thirds of the interstitial octahedral sites. Along the [0001] direction, consecutive O planes are separated by two Al 
layers. ZrC>2 displays three polymorphs at room temperature [19], either in its pure form (m-Zr02, monoclinic), 
or by dopage with substitutional cations (t-ZrC>2, tetragonal, and c-ZrC>2, cubic). The cubic structure corresponds 
to the CaF 2 lattice, which along the (111) direction leads to three possible terminations, labelled according to the 
composition of the two topmost layers: O-O- , O-Zr- and Zr-O-. t-ZrC>2 can be obtained from the cubic unit cell 
elongating one of the equivalent edges of the cube and introducing small opposite shifts at the O sublatticc positions. 
The m-Zr02 lattice is more complex, but along the [001] direction it defines a layered structure of chemical stacking 
sequence O-O-Zr. Here we have considered all low-index terminations of c-ZrC>2, and selected t-ZrC>2, m-Zr02, MgO 
and AI2O3 surfaces. Each surface is modelled as a periodic slab containing a vacuum region of 10 A. The width of 
the vacuum region is enough to inhibit interaction between neighbouring surfaces. The number of M^Oj, units in a 
slab depends on the particular structure, and ranges from 5 to 7. The slabs are symmetric about the central plane, to 
avoid an artificial divergent Madclung energy, although we have verified that all the conclusions of this work remain 
unaltered if asymmetric slabs are considered instead. Bulk-like behaviour is attained at the innermost central layers. 
The atomic positions are allowed to relax until the forces on the atoms are less than 0.06 eV/A. For comparison, a 
non-spin polarized (NSP) calculation is also performed for each case. The energy reduction due to spin polarization 
for all surfaces with a magnetic ground-state is shown in Table I. It has been obtained as the total energy difference 
between the relaxed SP and NSP slabs. The large values in the table are in the range of those obtained for magnetic 
bulk oxides [14], and should not be taken as an exact prediction of the actual values, but as a correct description of 
the order of magnitude and global trend. 

Table II summarizes our first main result. It shows the total magnetic moments at the surface O plane for all 
the structures which present a magnetic ground state, together with their in-planc (||) and normal (_L) p projections 
onto the surface. All spin moments arise almost exclusively from the p electrons at the valence band of the oxide. 
They reach large values, ranging from 0.8 to 1.6 [1b- Notice that only the O-rich oxide surfaces are magnetic. Both 
the cation-ended and the non-polar O surfaces not in the table do not show any spin polarization. The distinct 
characteristic of all these magnetic structures is that they are polar divergent surfaces [20] with a loss of cationic 
coordination for the outermost O atoms. This already suggests that the O spin polarization is intimately related 
to the decrease of the oxygen ionic charge at the surface with respect to the bulk. The correlation is evident after 
inspection of Table II, where the rightmost columns provide the O Mulliken populations at the surface layer and the 
bulk. In all these cases the surface O lose ionic charge approaching a charge neutral state. Moreover, the smaller the 
oxygen ionic charge, the larger its associated magnetic moment. 

Indeed, for the calculated surfaces with a non-magnetic ground state we obtain O charges close to the bulk ones. 
An example is given in table III, which compares the layer resolved Mulliken populations and moments for the two 
O-ended c-Zr02(lll) surfaces. ^From the results for c-Zr02(lll)o-o- i n table III, we can again observe the clear 
correlation between the reduction of the ionic charge and the size of the magnetic moment. This means that the 
spin polarization is a response of the system to the loss of transferred charge. Even if the magnetization is a local 
effect, rooted in the lack of donor electrons for the outermost O atoms, the interaction of these atoms with the lower 
adjacent planes cannot be neglected: the spin polarization at the surface layer induces a magnetic moment at the 
neighbouring planes, which decays as we deepen into the bulk. 

Our second main result is the finding that all magnetic surfaces are half-metallic. Half-metals are strong candidates 
for the development of spintronic devices, as they can naturally produce electronic currents of well-defined spin 
polarization [21]. Figure 1 shows the density of states (DOS) at the three uppermost layers and the bulk for two 
representative Al2O3(0001) and c-ZrO2(001) surfaces, both for SP (left) and NSP (right) calculations. In all the 
magnetic structures investigated, the Fermi level (E^?) crosses the valence band of the topmost layer, although some 
surfaces present a very localized minimum at Ep, as occurs in Al2O3(0001). Remarkably, the bands crossing F,p have 
always a well defined spin polarization, the majority spin states being full while the minority spin band is partially 
filled. The result, common to all the magnetic surfaces studied here, is a half-metallic system. 

The differences between the magnetic surfaces manifest in the details of the electronic distributions. This can be 
best observed regarding the shape of the DOS in figure 1. For c-ZrO2(001), the DOS at the topmost O plane is bulk- 
like except for the usual surface narrowing, both for the SP and NSP calculations. The effect of the spin exchange is 
mostly the shift of the minority spin band. At the SP Al2O3(0001) surface, the most significant effect is the reduction 
of the DOS at E^, something also observed at c-Zr02(lll)o-o- and MgO(lll). For these surfaces, the shape of the 
majority spin DOS is almost bulk-like, while a deep redistribution of charge affects the minority spin electrons. 
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The major changes observed in the minority spin DOS are also present in the spatial charge distributions. This 
can be directly seen in the left and central panels of figure 2, which depict the spin resolved charge density differences 
(CDD, total charge density minus the superposition of atomic charge densities) of the MgO(lll) and c-Zr0 2 (lll)o-o- 
surfaccs, within a plane normal to the surface along the bond direction between a surface O atom and a nearest 
neighbour at the layer below. The majority spin band being completely filled, as was shown in the DOS, the 
corresponding CDD adopts the spherical symmetry of the bulk. On the contrary, anisotropic lobular shapes appear 
in the minority charge. This assymetry between the spatial distributions of the spin states is a common feature to all 
the magnetic structures studied here, although the effect is specially relevant for the surfaces shown in figure 2. The 
rightmost panels of the figure show the corresponding spin density difference (SDD), CDD of majority spin minus 
CDD of minority spin. The crystal field induces a high directionality of the spin moments, involving orbitals along 
specific directions. For the c-Zr02(lll)o-O- an d MgO(lll) systems this leads to a highly anisotropic distribution 
of the magnetic charge, also evidenced in the m|| and m^ projections of table II. The specific dominant distribution 
depends on the particular structure: contained within the surface plane for c-Zr02(lll)o-o- and normal to the 
surface for MgO(lll). 

Finally, some considerations about energetics should be mentioned. To analyze the surface energetics requires 
to consider all mechanisms which stabilize the surfaces under study [4, 25], a huge task beyond the scope of this 
work. Nevertheless, it is well known that bulk metastable terminations coexist in grain boundaries, and they even 
become energetically favourable in ultrathin films and nanostructures [24]. Based on thermodynamic arguments 
[22, 23], the surface energy of the polar surfaces under study may be lowered below that of the corresponding non- 
polar terminations under high O pressures when the spin polarization is included. On the other hand, we have also 
performed preliminar studies of the effect of Zr vacancies in bulk Zr02 that confirm the existence of the magnetization 
mechanism presented here. Similar conclusions have been reported for bulk Hf02 in the presence of Hf vacancies [26] . 

In conclusion, we have proved the existence of local magnetism at the ideal O-ended polar surfaces of ceramic 
oxides. Large spin moments are formed at the outermost surface plane due to the lack of donor charge and the 
subsequent creation of 2p holes in the valence band of the oxide. This leads to a half-metallic ground state, with the 
majority spin band completely filled and only minority spin states at the Fermi level. The orientation of the magnetic 
moment with respect to the surface plane can be tuned selecting adequate crystal orientations. The fact that polar 
surfaces of nonmagnetic materials with large bulk bandgap exhibit a 100 % spin polarization at the Fermi level opens 
a new route to the manipulation of spin currents. In fact, the mechanism presented here can explain the unexpected 
magnetism measured in thin films of Hf02 and Zr02 [8, 9], relating it to the existence of cation defects. Even more, 
the present ability to produce ultrathin films and nanometric grains of these materials make feasible the fine control 
of their structure and stoichiometry, implying the possibility to extend this surface magnetism to any low-dimensional 
configuration. 
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TABLE I: Energy reduction (eV /M x O y unit) of the SP relaxed slab with respect to the corresponding NSP case for all magnetic 
surfaces under study. 
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TABLE II: Spin moments (in /is) at the topmost layer of the surfaces of Table I, together with the decomposition of the 
p orbital contribution parallel and normal to the surface plane. The rightmost columns provide the corresponding Mulliken 
charge populations (Q) compared to the inner bulk value (Qf,). The Oi and O2 entries for the non-cubic Zr0 2 structures refer 
to the two inequivalent in-plane positions. 
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FIG. 1: DOS for the three topmost planes and the bulk of Al 2 O 3 (0001) and c-ZrO 2 (001) after SP (left) and NSP (right) 
calculations. Positive (negative) values correspond to majority (minority) spin states. Energies are referred to the Fermi level. 
The total spin moments at each layer are indicated for the SP case. 
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TABLE III: Layer resolved spin moments (in /is) and Mulliken populations for the O-ended c-Zr02(lll) surfaces. The layers, 
labelled according to their composition, are ordered from the surface (top row) to the bulk. 
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FIG. 2: Left and middle panels: spin resolved CDD for the (top) MgO(lll) and (bottom) c-Zr02(lll)o-o- surfaces. Right 
panels: corresponding SDD. 
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